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ABSTRACT 

We present the long-term monitoring of the High Mass X-ray Binary GX 301-2 performed with the Super- 
AGILE instrument on-board the AGILE mission. The source was monitored in the 20-60 keV energy band 
during the first year of the mission from 2007 July 17 to 2008 August 31, covering about one whole orbital pe- 
riod and three more pre-periastron passages for a total net observation time of about 3.7 Ms. The SuperAGILE 
dataset represents one of the most continuous and complete monitoring at hard X-ray energies of the 41.5 day 
long binary period available to date. 

The source behavior was characterized at all orbital phases in terms of hard X-ray flux, spectral hardness, spin 
period history, pulsed fraction and pulse shape profile. We also complemented the SuperAGILE observations 
with the soft X-ray data of the RossiXTE/ASM. Our analysis shows a clear orbital modulation of the spectral 
hardness, with peaks in correspondence with the pre-periastron flare and near phase 0.25. The hardness peaks 
we found could be related with the wind-plus-stream accretion model proposed in order to explain the orbital 
light curve modulation of GX 301-2. 

Timing analysis of the pulsar spin period shows that the secular trend of the ~680 s pulse period is consistent 
with the previous observations, although there is evidence of a slight decrease in the spin-down rate. The 
analysis of the hard X-ray pulsed emission also showed a variable pulse shape profile as a function of the 
orbital phase, with substructures detected near the passage at the periastron, and a clear modulation of the 
pulsed fraction, which appears in turn strongly anti-correlated with the source intensity. 
Subject headings: pulsars: individual: GX 3012 stars: neutron X-rays: binaries 



1. INTRODUCTION 

GX 301-2 is a High Mass X-ray Binary (HMXB) sys- 
tem containing an X-ray pulsar, with a ~685 s rotation pe- 
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riod, orbiting the B-emission line hypergiant star (Bl Ia+) 
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with a balloon observation, subsequently 
discovered the 700 s X-ray pulsation. The 
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a period of -41 . 5 days dWhite et al.lll978t ISato et alj 119861 : 
iKohet all \l99% iDoroshenko et al.l 120081) The mass func- 
tion of the system is 31.1 M Q with a companion mass in the 
range 39 < M < 53 M , estima ted by the measur ement of 
optical radial velocity amplitud e dKaper et alj|2006l) . The ra- 
dius of Wray 977 obtained by iKaper et alj ( 120061) by fitting 
atmosphere model is 62 R G , while the effective temperature is 
about 1.8 x 10 4 K. 

GX 301-2 shows regular X-ray flares about two days before 
the periastron passage, while a secondary inten sity peak was 
some t imes observed near the apastron passage dPravdo et alj 
[1991 iKoh et alJ[l997l iPravdo & GhosMl200lh. A ssuming a 
distance of 3 kpc, as derived bv lKaper et alj d2006l) . the source 
luminosity during pre-periastron (PP) flares reaches values 
of ~ 10 37 erg s , about 25 times larger than the luminosity 
showed by the source in the other orbital phases. 

Since its discovery, GX 301-2 exhibited unpredictable 
torque variations superimposed on a secular spin trend, 
which passed from a spin-up to a spin-down state between 
MJD 48500 and MJD 49500. Two rapid spin-u p episodes 
with v ~ 2-5 x 10~ 12 Hzs~' were also found by IKoh et alj 
(119971) analyzing BATSE data, th us confirm i ng the erratic be- 
havior of the pulsar spin period. IKoh et all (11997b suggested 
that the long-term spin-up trend observed between 1984 and 
1994 might not be regular and continuous, but due to sev- 
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eral rapid spin-up episodes, which in turn c ould be related to 
the for mation of a temporary accretion disk. iPravdo & Ghoshl 
(12001 . using BATSE data, found evidence that the hard X-ray 
(20-50 keV) intensity of the binary system during quiescent 
orbital phases is correlated with the pulse period, suggesting 
a possible relation between the spin period behavior and the 
mass accretion rate. 

Several wind accretion models were proposed to explain 
both the regular PP and apastron flaring activities of the 
source. A dense disk around Wray 977, which intercepts the 
neutron star (NS) orbit near the peria stron and the apastron, 
was p roposed by IPravdo et al.l (11995b and IPravdo & Ghoshl 
( 12001b . However, this model seems unable to provide a 
satisfactory modeli ng of the o bserved orbital X-ray light 
curve of the source dLeahy||2002l) and optical observations by 
iKaper etafl (12006b found no confirmation of the presence of 
a disc. The existence of a high density stre am of matter i n 
addition to the stel l ar wind was proposed b y lHaberll (11991b : 
ILei^(fl99ll 12002): Leahy & Kostka (2008). The stream out- 
flows from Wray 977 due to tidal interaction intercepts the NS 
orbit near the periastron and the apastron, causing a large in- 
crease in the mass accretion rate with a sudden enhancement 
of the X-ray luminosity from the neutron star. Signatures of 
the pr esence of such a stream were found by iSaraswat et alj 
d 1996b analyzing the low-ener gy excess in the AS CA spec- 
trum of GX 301-2, and also bv lKaper et ail d2006b from the 
orbital modulation of the spectral lines detected in the stellar 
wind of the optical companion. The wind-plus-stream model 
predicts also a large variation in the column density as a func- 
tion of the orbital phase. The Nn is expected to change from a 
value of ~ 10 23 cirT 2 at orbital phase 0.8 to values greater than 
10 24 crrT 2 during the PP flares and between orbital phases 
0.2-0.3. 

We observed the HMXB GX 301-2 with the SuperAGILE 
(SA) experiment onboard the Astro-rivelatore Gamma ad Im- 
magini LEggero (AGILE) satellite during the first year of the 
mission. AGILE, launched on 2007 April 23, is a scientific 
mission of th e Italian Space Age ncy dedicated to high-energy 
astrophysics dTavani et al.l2009h . In the following sections we 
report on the analysis of the SuperAGILE data, in the energy 
range 20-60 keV. In Section [2] we describe the SuperAGILE 
experiment and the observations of GX 301-2 we performed, 
in Section[3]we present the timing analysis of the pulsar emis- 
sion and in Section |4] we analyze the orbital light curve of 
GX 301-2 using the hard X-ray SuperAGILE and the soft X- 
ray RXTE/ ASM data and the orbital variability of the source 
pulsed fraction. Discussion of the data analysis and the in- 
terpretation of the results are presented in Section|5] while in 
Section|6]we draw our conclusions. 

2. SUPERAGILE OBSERVATIONS AND DATA REDUCTION 

2.1. The SuperAGILE Experiment 

SuperAGILE (see iFeroci et"afl 120071 for an extensive de- 
scription of the experiment) is the hard X-ray monitor of the 
AGILE mission and is based on four one-dimensional coded- 
mask instruments, each one equipped with a ~ 360 cm 2 sili- 
con microstrip detector and a tungsten mask-collimator sys- 
tem. The experiment has ~1 sr field of view (FoV), 6 arcmin 
angular resolution, and a point source location accuracy bet- 
ter than 2 arcmin for intense sources. SuperAGILE is able 
to detect X-ray photons between 18 keV and 60 keV with a 
resolution of ~8 keV FWHM and an on-axis effective area of 
250 cm 2 . The time resolution of the instrument is 2 ps while 



the accuracy is about 5 ps. The dead time is negligible in nor- 
mal operation conditions, except for the passages through the 
South Atlantic Geomagnetic Anomaly where the vetoes from 
the AGILE anticoincidence systems reduce the SA live time 
to ~10% or less. In the standard analysis this fraction of the 
orbit is usually removed from the data. 

In standard operation mode SuperAGILE provides event by 
event data, from which position and normalized count rates 
for each source detected in the FoV are derived. The data are 
downloaded from the spacecraft about every 100 minutes dur- 
ing the AGILE passage over the Malindi ground station, sent 
to the ASI Science Data Center (ASDC) in Rome and then 
processed by the SuperAGILE Scientific Pipeline (SASOA) 
developed by the Supe rAGILE team at INAF/IASF Rome 
dLazzarotto et al. 1 120081) and running both at the ASDC and 
at IASF. The automated analysis of the SuperAGILE data au- 
tonomously updates a source catalog and orbital light curves 
for bright sources are then posted at the ASDC web sit^ZJ 

2.2. SuperAGILE Observations ofGX 301-2 

SuperAGILE observations of the HMXB GX 301-2 started 
in mid 2007 July, just after the end of the AGILE Commission- 
ing Phase, as part of the Science Verification Phase. Indeed, 
GX 301-2 was the first light of the experiment. Due to the so- 
lar constraints and the AGILE pointing strateg)0 the satellite 
is operated by performing long observations, typically of 2- 
4 weeks duration, during which the pointing direction slowly 
drifts (at a rate of ~1° day -1 ). SuperAGILE performed six 
long observations of GX 301-2, for a total net exposure time 
of about 3.7 x 10 6 s. Table Q] shows the journal of the obser- 
vations. The first observation and part of the second one were 
performed during the AGILE Science Verification Phase that 
started in mid 2007 July and ended on 2007 November 30. 
The other observations were carried out during the AGILE 
AOl program, which started on 2007 December 1 and lasted 
one year. The last column in Table Q] shows the source off- 
axis angle in the SuperAGILE FoV for each observation. The 
configuration of SuperAGILE implies that the source counts 
collected by t he instrument stro ngly depend on the off-axis 
angle (see e.g. lFeroci et al.ll2~007b . and thus also the detection 
significance is dependent on the source position in the FoV. 

The SuperAGILE observations took place in different seg- 
ments of the GX 301-2 binary orbit. In Figure Q] we graph- 
ically show the orbital phases of our six observations: three 
of them covered a large fraction of a complete binary orbit 
(Observations ID II, III and V), while the other three were 
shorter, and covered the equivalent of one full binary orbit all 
together, but arranged over three different orbits. The Super- 
AGILE data thus allow for a often continuous monitoring of 
the complete 41.5 days binary cycle. It is worth noting that 
the monitoring of four different passages at the pre-periastron 
orbital phase (0.95) was included in the SuperAGILE obser- 
vations. 

2.3. Data Reduction and Analysis 

SuperAGILE data r eduction was performed using the stan- 
dard SASOA pipeline dLazzarotto et al.l l2"008l). which extracts 
high level products from the photon-by-photon SuperAG- 
ILE data. Production of source images was performed with 
the SuperAGILE Enhanced Multi Imaging (EMI) procedure 

21 http://agile.asdc.asi.it/sagilecat_sources.html 

22 http ://agile. asdc . asi. it/current. pointing. html 
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Fig. 1, — Graphical representations of the GX 301-2 binary system. The solid thick line in each panel shows the orbital phases observed by SuperAGILE during 
the six observation blocks summarized in TablefT] Ticks on the Neutron Star orbits represent the orbital phases with step 0.1 as displayed in the upper-left panel, 
the NS orbital motion is counter-clockwise. Units along the axis are in solar radii while the arrow in the upper-left panel represents the observer direction. Roman 
numerals in each panel indicate the Observation ID. 



TABLE 1 

Journal of the AGILE Observations of GX 301-2. 



Obs. ID 


Mission Phase 


MJD Start 


MJD Stop 


Orbital Phase 


Off-axis Angle (deg) 


I 


SVP 


54298.4 


54305.3 


0.49-0.66 


7-10 


II 


SVP, AOl 


54311.6 


54339.1 


0.81-0.47 


4-15 


III 


AOl 


54473.6 


54497.7 


0.71-0.29 


5-18 


IV 


AOl 


54510.6 


54527.3 


0.60-0.01 


10-17 


V 


AOl 


54647.6 


54672.9 


0.91-0.52 


14-20 


VI 


AOl 


54694.3 


54709.5 


0.03-0.40 


12-17 



Note. — First column indicates the Observation ID, the second column specifies the AGILE mission phase: SVP=Science Verification Phase, AOl =Scientific 
Observations Program pointings (see AGILE Mission Announcement of Opportunity Cycle- 1 : http://agile.asdc.asi.it/). 



dEvangelista et a l. 2008) which allows fast and accurate anal- 
ysis for bright sources. For each AGILE orbit the source 
count rate was automatically extracted in the energy range 20- 
60 keV, excluding the time spent by the satellite in crossing 
the South Atlantic Geomagnetic Anomaly and the time when 
the Earth occulted the line of sight to the source. Normalized 
count rates (in units of counts cirT 2 s~') were then obtained by 
considering the effective area of the exposed po rtion of the de- 
tector s and assuming a Crab-like spectrum (e.g. lFrontera et al] 
120071) . Fluxes in unit of the Crab flux were calculated by ap- 
plying a normalization factor of 0.15 to the normalized count 
rates. Uncertainties on the flux measurements were estimated 
by adding to the statistical e rror a systematic e rror equal to 
the 8% of the measured flux dFeroci et a l. 2009). In order to 
avoid the introduction of biases in the light curve we did not 
select the source detections for significance but all the obser- 
vations were used in the analysis. A detailed description of 
the method used for t he generation of the SuperAGILE light 
curves can be found in lFeroci et al.l (120091) . The resulting 20- 
60 keV light curve of GX 301-2 with a ~ 6 x 10 3 s binning is 
shown in Figure [2] 

SuperAGILE events are time-tagged with microsecond pre- 
cision by the AGILE on-board Pulse-Per-Second (PPS) sys- 
tem, which is synchronized to UTC once a second by the 



Global Positioning System (GPS). In order to perform timing 
analysis all the SuperAGILE events are converted on ground 
to the Terrestrial Dynamical Time (TDT) reference system. 
The photon arrival time (TOA) is then converted to Bary center 
Dynamical Time and the correction for the arrival delay at the 
Solar System Barycenter is calculated by using the JPL plan- 
etary ephemeris DE405. A more de tailed description of th e 
correction algorithm may be found in lPellizzoni et al.l d2009l) . 
Correction for the time delay introduced by the binary orbital 
motion was performed using the Blandford-Teukolsky m odel 
dBlandford & Teukols]gvl[T976l iTavlor & Weisberd[l989l) . To 
perform the TOA correction w e used the orbital parameters 
provided by iKoh et all ( 119971) . except for the orbital period 
(P 01 b) and the time of periastron passag e (To), for which 
we uti lized the updated values derived by iDoroshenko et al.l 
d2008l) . These parameters are given in Table [2] The quoted 
uncertainties are given at 1 cr level throughout the paper, un- 
less otherwise specified. 

3. ANALYSIS OF THE SPIN PERIOD 

3.1. Pulse-shape Profile 

A search for the pulsed emission from the NS was per- 
formed on the corrected events of all the SuperAGILE source 
observations using the standard epoch folding technique 
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Fig. 2. — Light curve of GX 301-2 as observed by SuperA GILE. The source was mo nitored for MJD 54298 to MJD 54709. Dashed lines indicate the periastron 
passages on the basis of the orbital parameters calculated by[Doroshenko et al. 1 2008ft . 



TABLE 2 
Orbital Parameters of GX 301-2. 



Parameter Koh et al. ( 1 997) 
P orb (day) 41.498 ±0.006 
e 0.462 ± 0.042 

368.3 ± 11.1 

310.4 ±4.2 
48802.79 ± 0.36 



Doroshenko et al. (2008) 



a r sin i (lt-s) 
oj (deg) 
To (MJD) 



41.482 ±0.006 



53531.63 ±0.03 



Note. — Param eters are derived by IKoh et all 119971) and by 
Doroshenko et al. 1 2008). The uncertainties are at the 3cr level both for the 
Koh and the Doroshenko parameters. 



dLeahv et al.lll983l) . Due to the different off-axis angles to 
which the source was observed in the different pointings (see 
Table [TJ, corresponding to significantly different exposed ef- 
fective areas, for each orbital phase bin we selected the avail- 
able data with the highest signal-to-noise ratio, in order to 
obtain the best determination of the pulse shape. In case of 
pulse shape variability between different orbital cycles, this 
choice biases the pulse shape to the specific selection. How- 
ever, we verified that in the other SuperAGILE observations 
of the same orbital phase bin the pulse shape was consis- 
tent, within the statistical uncertainties. Following this ap- 
proach, pulse shape profiles for orbital phases between 0.26 
and 0.64 and between 0.81 and 0.06 were obtained from data 
collected during Observation I and II, while the other pulse 
shape profiles were built using data taken during Observation 
III. The width of the phase intervals was chosen according to 
the signal-to-noise ratio in the profiles. This implies to have 
a phase bin equal to ~0.10 (~4 days) for the orbital phase 
interval 0.06-0.34, where the source flux is lower, while a 
phase bin size <0.05 was used for the other orbital phases 
with higher signal-to-noise. We were not able to detect NS 
pulsation in orbital intervals 0.16-0.26, 0.34-0.38, 0.43-0.49, 
0.54-0.59 and 0.64-0.71. 

Due to the erratic behavior of the spin period of th e wind- 
fed pu l sar even on sh o rt tim e scales, as reported bv | Nagasd 
(11 9891) . iBildsten et all (119971) . iPravdo & Glioshl (120011), and 
taking into account the errors on period determination, we 
were not able to determine the absolute phasing of the pulse 
shape profile. Hence, we cross-correlated each pulse pro- 
file with a template obtained in the orbital phase 0.96-0.01 



(lower-right panel in Figure [3), which represents the pulse 
shape with the highest signal-to-noise. The resulting phase 
shifts were then used to build the background subtracted 
phase-aligned folded light curves shown in Figure [3] De- 
spite the different statistical quality of the curves, a variability 
along the orbital cycle can be clearly identified. The main 
peaks change in duration, shape and height. The main peak, 
at spin phase 0.0, shows a significant double peak structure 
only at orbital phases between 0.96 and 0.06, where the most 
dramatic flux evolution of the source occurs. Although these 
are the orbital phases where the statistical quality is higher, 
the insurgence of the minor peak does not appear as an obser- 
vational bias: the evolution is gradual and continuous in the 
three panels between phases 0.91 and 0.06, all with compara- 
ble quality. 

Another variability feature that is worth remarking is the 
behavior of the interpulse regions. At spin phase around 0.7 a 
small peak appears at orbital phase 0.88-0.91. Its onset seems 
already anticipated in the previous orbital phase intervals. The 
other interpulse region, around spin phase 0.3, shows instead 
an intensity variability. The source flux in this interval is sig- 
nificantly higher, with respect to the peaks, in the spin phase 
interval between 0.96 and 0.06, the same where the small peak 
mentioned above appears. 

3.2. Secular Evolution of the Spin Period 

The SuperAGILE long-term and nearly continuous moni- 
toring of GX 301-2 allowed us to measure the spin period 
at different epochs and put it in the context of the secu- 
lar timing evolution of the source, as available from the lit- 
erature. Figure Ufa) shows the spin period P as a func- 
tion of the time from the discovery of the pulsed emis- 
sion by IWhite et all d 19761) to the last SuperAGILE obser- 
vation (MJD 54709.5). This plot confirms that GX 301-2 
is subject to unpredictable torque variations superimposed 
on long-term smooth variations of the spin period. Dur- 
ing the last ~30 years the spin period changed passing from 
a spin-up state, during which P decreased from ~700.5 s 
to ~675.4 s, to a spin-down condition reaching the current 
value of P ~687.3 s. The Figure clearly shows th e two 
rapid spin-up episodes identified by IKoh et all ( 11997b ana- 
lyzing CGRO/BKYSE observations between MJD 48440 and 
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Fig. 3. — GX 301-2 pulse shape profiles. Each panel represents the pulse shape profile, after the background subtraction, in an orbital phase interval as reported 
on top of each figure. See the text for a detailed description of the methods used to obtain the pulse shape profiles. 



MJD 48463 and between MJ D 49230 and MJD 49245. C om- 
paring BA TSEp^l RXTEfPCA jKrevkenbohm et al.l l2004l). IN- 
TEGRAL dDoroshenko et alj|2008l) and SuperAGILE data in 
Figure Ufa) we find evidence of at least two other spin-up 
episodes, occurred between the observations at MJD 51830 
and MJD 52789, and around MJD 54000. 

Figure HJb) displays the secular variation of the GX 301-2 
pulse period during the SuperAGILE observations. P varies 
in the range 685.3-688.0 s in about 410 days. We find a pe- 
riod evolution consistent with the constant spin down rate of 
4.0 + 0.3 x 10 -8 s s _1 . However, the short term variability of 
the period clearly shows significant timing noise around the 
general trend, and the linear fit (dashed line in Figure H|b)) is 
indeed a poor modeling of the period evolution, leading to a 
reduced^ 2 value of 2.9 (34 dof). 

23 http ://gammaray. msfc . nasa. gov/batse/pulsar 



4. ORBITAL VARIABILITY 

4.1. SuperAGILE and RXTE/ASM Folded Light Curves 

Using all the SuperAGILE data collected during the obser- 
vations listed in Table Q] we obtained the complete, folded 
orbital light curve of GX 301-2 shown in the top panel of Fig- 
ure [5] The orbital parameters (P or b, To) used in the folding 
procedure were the same used for the correction for the time 
delay discussed i n Section [3] The differe nce in the orbital 
solutio n found bv lDoroshenko et all d2008l) and bv lKoh et alj 
(1 1997b results in a 5% phase shift of the PP flare in the or- 
bital folding, while the shape of the folded light curve does 
not change noticeably. In order to obtain a folded light curve 
in the soft X-ray band we used public data f rom the All 
Sky Monitor onboard RXTE dLevine et al.lll996l) . ASM data 
for GX 301-2 cover the time span from MJD 50087.285 to 
MJD 54720.745 and consist of source count rates for each 
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Fig. 4. — (a) Secular variations of GX 301-2 p ulse period. Dat a 
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iDoroshenko et al. 2008), filled squares in highlighted region (SuperAGILE 
data), (b) Zoomed view of the SuperAGILE data. The dashed line represents 
the best-fit secular variation of the pulse period assuming a linear trend. 



ASM dwell (1 dwell = 90 s observation) in the band 2 keV- 
12 keV. The dwell -by-dwell data were filtered according to 
the criteria reported on RXTE/ ASM web sitePI and then a light 
curve was built from the cleaned data. The ASM folded light 
curve for GX 301-2 is displayed in the second panel of Fig- 
ure |5](solid line). 

Both the SuperAGILE and the ASM light curves are 
given in units of mCrab. We used the standard value of 
1 Crab = 75 counts s as a conversion factor for the 
ASM count rates, while the SuperAGILE normalization is 
0.15 counts crrT 2 s _1 . The two folded orbital light curves 
clearly show the PP flaring activity, with an intensity peak 
observed at orbital phase between ~0.9 and ~0.0 and a full 
width at half maximum of about 0. 1 in phase. Fitting the PP 
flare by using a Gaussian function plus a constant term, we 
obtained a value of orbital phase of 0.954±0.005 for the flare 
position, with a FWHM value of 0.109. The remaining or- 

24 http://heasarc.gsfc.nasa.gov/docs/xte/SOF/asmlc.html 
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Fig. 5. — Top panel: SuperAGILE orbital light curve (20-60 keV) of 
GX 301-2 folded at orbital period P orb =41.482 d. Second panel: RXTE/ ASM 
orbital folded light curve (2-12 keV) of GX 301-2 obtained using the whole 
ASM dataset from MJD 50087.285 to MJD 54720.745 (black solid line) and 
only the observations contemporaneous to SuperAGILE (red stars). Both the 
folded light curves are given in mCrab units. Third panel: Hardness ratio 
(20-60 keV to 2-12 keV) calculated using the whole SuperAGILE and ASM 
dataset (filled circle) and the contemporaneous observations only (red stars). 
Bottom panel: pulsed fraction (Pf rac ) as a function of the orbital phase. Thick 
points represent the mean pulsed fraction values in the orbital phases 0.05- 
0.39, 0.39-0.86 and 0.86-0.05. The sharp P ft ac minimum which corresponds 
with the pre-periastron flare is clearly visible. 

bital interval exhibits a different behavior in the hard and soft 
X-ray data: in the 20-60 keV energy band a nearly constant 
persistent emission is detected between phases 0.1 and 0.8, 
while in the same phase interval the soft X-ray data show a 
broad hump, covering phases between 0.5 and 1.0. 

4.2. Variability of the pulsed fraction 

The variability in the shape of the pulse as a function of 
the orbital phase shown in Figure [3] motivated us to study the 
pulsed fraction across the binary orbital cycle, to determine its 
possible variability and correlation with other source param- 
eters, such as flux and hardness. We operationally define the 
pulsed fraction (Pf ra c) as the ratio between the pulsed counts, 
as derived directly from folding the event list, including pho- 
tons from the background and from other sources in the FoV, 
and the total counts determined by the imaging procedure, in- 
cluding the counts from GX 301-2 only. We performed the 
analysis of the pulsed fraction as a function of the pulsar or- 
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TABLE 3 

Mean Hardness Ratio of the PP Flares 



MJD Start 


MJD Stop 


Phase Start 


Phase Stop 


HR 


54315.84 


54319.37 


0.91 


0.99 


10.4 ± 0.6 


54481.77 


54485.30 


0.91 


0.99 


17.5 ± 1.2 


54523.25 


54526.78 


0.91 


0.99 


21.8 ± 1.7 


54647.70 


54652.22 


0.91 


0.99 


14.0 ± 0.9 



bital phase for all the available SuperAGILE observations dur- 
ing which the pulsation was detected. The results are shown 
in the bottom panel of Figure [5] where Pf lac is plotted as a 
function of the orbital phase. 

ff rac exhibits large variations across the binary cycle, pass- 
ing from a minimum of Pf rac = 0.15 + 0.02 at orbital phase 
~0.0 to a maximum of Pf rac = 0.64 + 0.12 at orbital phase 
~0.5. With the aim of quantifying the significance of the vari- 
ability of the pulsed fraction, we took the average of ff rac in 
the phase intervals where it appears consistent with a con- 
stant value, that is 0.05-0.39, 0.39-0.86 and 0.86-0.05. Thick 
red points in the bottom panel of Figure |5]represent the mean 
Pf rac values calculated during the three selected orbital inter- 
vals. The values in the three phase intervals are 0.188+0.009, 
0.387+0.025 and 0.556+0.034, ordered by phase. This analy- 
sis shows that it is possible to identify three different "states" 
in the pulsed fraction as a function of the NS orbital motion, 
with the intermediate one occurring just after the PP flares 
and lasting approximately up to the phase corresponding to 
the superior conjunction of the binary system. 

The bottom panel in Figure [5] shows that the pulsed frac- 
tion reaches its minimum at the phase when the large flare 
occurs, that is where the largest variations in the pulse shape 
are detected. We then investigated the possible dependence 
of ff m c on the source flux. The results are shown in Fig- 
ure |6ja), where the value of the pulsed fraction against the 
source flux is displayed, independently of the orbital or spin 
phase of the source. A clear anti-correlation between these 
two quantities is found, confirming it as a general feature 
and not only an event occurring during the PP periodic flare. 
The plot shows that the pulsed fraction decreases rapidly at 
increasing fluxes, reaching a plateau for large fluxes corre- 
sponding to a pulsed fraction of 15%, approximately a factor 
4 smaller than the maxi mum value. A similar behavior was 
also recently reported bv lLutovinov & Tsygankovl (120091) us- 
ing INTEGRAL/IBIS data. 

4.3. Hardness Ratio Versus Orbital Phase 

Given the different orbital flux modulation at soft and hard 
X-rays (Figure |5), we investigated the period-averaged time 
behavior of the source hardness along the orbital light curve. 
To this aim we computed an hardness ratio (HR) by using the 
flux measurements by SuperAGILE and ASM, that is between 
the energy ranges 20-60 keV and 2-12 keV. The third panel 
in Figure[5]shows the result. Filled circles indicate the orbital 
phase-resolved hardness ratio obtained by using the complete 
temporal data set from both experiments. Moreover, star-like 
symbols provide the same quantity, but computed after select- 
ing the ASM time intervals coeval with the SuperAGILE data 
sets (the ASM light curve with this selection is shown using 
the same symbols in the second panel of the same figure). 

The time profile of the HR computed using both time se- 
lections indicates a clear modulation of the source spectral 
hardness as a function of the orbital phase, where we iden- 
tify the following features. A spectral hardening is detected 
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Fig. 6. — (a) Pf rac as a function of the source intensity. The clear anti- 
correlation between pulsed fraction and source flux is ascribed to the different 
amount of matter which accumulates near and accretes onto the neutron star, 
(b) and (c) Hardness ratio as a function of the source flux in the 20-60 keV 
(b) and 2-12 keV (c) energy bands. 

during the pre-periastron flare: the hard X-ray flux rises faster 
than the soft X-rays, and lasts longer. Then, when the soft 
X-ray flux drops at its minimum orbital level, around phase 
0.2, the hard X-rays hold their intensity level, indicating an 
increasing hardness of the emitted radiation. Subsequently, 
the rise of the soft X-ray flux after phase 0.5 (the apastron) 
is not accompanied by a corresponding increase at hard X- 
rays, indicating a softer emission in this part of the orbit. We 
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note that using the simultaneous data only (stars in the HR 
plot), the above features are confirmed and, within the lim- 
its of the lower statistical quality, emphasized. Actually, the 
source spectral hardness changes also from one binary orbit 
to another. In Table [3] we provide the HR values computed in 
the orbital phases including the PP flare in the four available 
SuperAGILE observations. The value of the HR changes up 
to a factor of 2. 

5. DISCUSSION 

The pointing strategy of the AGILE mission offered Su- 
perAGILE the opportunity of studying the temporal proper- 
ties of GX 301-2 in the hard X-rays over repeated ~month- 
long stretches of time, covering six 41.5 days binary orbits 
at different phases. We observed a modulation of the hard 
X-ray flux along the orbit in general agreement with what 
reported over the past two decades by previous experiments 
in soft and hard X-rays. Fitting the source spectrum in the 
20-60 keV energy band with a cut-off power-law with pho- 
ton index fixed to F=1.05 ( La Barbera et alj|2005l) . we obtain 
£ cut =24.4+|-2 keV and E fold =7.8+<^ keV for the PP flare oc- 
curred around MJD 54316. Usin g these parameter s, and as- 
suming a source distance of 3 kpc dKaper et alJ2006l) . we con- 
vert the SuperAGILE counts to physical units obtaining a lu- 
minosity of ~2.3xl0 37 erg s _1 . Under the same assumptions, 
the luminosity during the quiescent orbital phases is approxi- 
mately one order of magnitude smaller. 

The only complete folded orbital light curve of GX 301- 
2 at energies above 20 keV available in the literature, to 
which we can compare the SuperAGILE results, was built 
from the data of the BATSE experiment , using the Earth 
occul tation technique dPravdo et al.l 1 19951 IPravdo & Ghoshl 
12001b or the pulsed flux measurements dKoh et al.l 119971) . 
The SuperAGILE data show significant differences. In con- 
trast to the BATSE ob servations in the same energy range 
(IPravdo & Ghoshll2001l) . the SuperAGILE data do not show 
any secondary peak in the folded light curve near the passage 
at the apastron (orbital phase ~0.5). This second peak was in- 
deed weakly detected in the BATSE m easurements of the total 
source flux (IPravdo & Ghosh ] |2001h. while i t was outstand- 
ing in the pulsed flux only (IKoh et all 1 19971) . The different 
shape of the orbital light curve of the total flux as measured 
by BATSE and SuperAGILE may be ascribed to intrinsic time 
variability between the different orbits and/or epochs of the 
observations (1991-1997 for BATSE and 2007-2008 for Su- 
perAGILE) and/or to the different effective bandpasses of the 
two experiments (harder for BATSE than for SuperAGILE de- 
spite their similar nominal energy range). Actually, hints to- 
ward time variability near the apastron phase were already 
provided b y the epoch-resolved BA TSE folded light curves 
reported by IPravdo & Ghoshl ( 120011) . where only two out of 
four curves show a broad hump around phase 0.5, although 
the statistical quality of the data does not allow for a definite 
statement. Also the recent Swift/BAT da ta of three orbital cy- 
cles reported by iPestalozzi et alj J2009) show bumps around 
phase ~0.5 in two contiguous orbital cycles around 2008 De- 
cember (unfortunately these time periods were not covered by 
SA observations). 

Instea d, the pronounced bump around the apastron phase 
found by lKoh et al.l (119971) in the pulsed flux appears consis- 
tent with our results, when the behavior of the pulsed frac- 
tion along the binary orbit is considered. As we show in Fig- 
ure|5] we found a large orbital variation of the pulsed fraction: 
within the statistical accuracy of our data, the pulsed fraction 



exhibits a step-like behavior along the orbit. The minimum 
value is reached at the periastron, while a broad maximum, 
nearly 3 times as high, is reached at phases between 0.4 and 
0.9. This implies that the pulsed flux increases in the phase 
range corresponding to the apastron, coherently with what 
was found in the BATSE pulsed data. 

Interestingly, the same orbital phases where the pulsed frac- 
tion is highest are also those where the emitted radiation is 
softest. In fact, we analyzed the behavior of the source hard- 
ness along the orbit by using the SuperAGILE data in com- 
bination with the public RXTE/ ASM data. We found that the 
radiation emitted by the source is harder near the periastron 
and at phases around 0.2-0.3, while it becomes significantly 
softer at phases 0.5-0.8 (Figure [5]). The hardening around 
phase 0.2-0.3 was also detected by IPravdo & Ghoshl ( 120011) 
through an analysis similar to ours, using the BATSE and 
ASM data, but they did not report evidence for the harden- 
ing near the periastron, that is instead the main feature in our 
data. More in general, our data suggest an anti-correlation be- 
tween the hardness and the pulsed fraction at hard X-rays (as 
well as between the pulsed fraction and the flux, Figure|6ta)). 

The general behavior of the source is summarized in the 
synoptic panels of Figure [7j where the parameters are graphi- 
cally displayed along the binary orbit. These plots show that 
even if the only significant feature in the flux modulation ob- 
served by SuperAGILE occurs at the pre-periastron phase, 
the emission mechanism undergoes important changes also 
in other parts of the orbit. The passage at the periastron is 
accompanied with a hardening of the radiation and with a 
sudden and large decrease of the p ulsed fraction, as already 
reported by La Bar bera et al.l (120051). using the B eppoSAX ob- 
servations, and iLutovinov &Tsygankovl(l2009l) using INTE- 
GRAL data. The pulsed fraction shows (Figure [6| a)) a clear 
anti-correlation with the source flux, thus presumably with 
the accretion rate. Instead, the hardness ratio does not show 
a clear dependence on the X-ray flux (and then the accretion 
rate), suggesting a more complex relation. It is interesting 
to note that the largest drop in the pulsed fraction is approx- 
imately a factor of 3, from ~0.55 to ~0.2, while the flux in- 
crease in the same phase interval is about a factor of 10. This 
implies that the excess hard radiation emitted during the pre- 
periastron flare does not represent a proportional luminosity 
increase of the same regions emitting in the preceding orbital 
phases, that would leave the pulsed fraction unchanged, nor 
an unrelated emission from a site not subject to occultation by 
the NS upon spinning, that would cause a larger drop in the 
pulsed fraction. Thus, either the site of the excess emission is 
only partially occulted to the line of sight as the NS rotates, 
or it is associated with an intervening radiative transfer mech- 
anism able to sm ear the coherence of the photon arrival times 
(as proposed by iMakino et al.lll985l) . Moving from the pe- 
riastron to the apastron, the pulsed fraction seems to recover 
gradually to its maximum value. Instead, the hardness ratio 
shows a further, smaller peak at phases near 0.25 mostly due 
to a depression in the soft X-ray flux. Thus the cause inter- 
vening in obscuring or depressing the emission of soft X-rays 
where the NS approaches the apastron phase, leaves the per- 
sistent hard X-ray component unaffected, as witnessed by the 
u nchanging flux level mea sured by SuperAGILE. 

iLeahvl ( 120021) (see also iLeahv & Kostkall2008h quite con- 
vincingly modeled the orbital modulation of the soft X-ray 
emission in the GX 301-2 system as driven by wind accre- 
tion plus a mass stream, drawing an Archimedes spiral path 
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(c) 

Fig. 7. — Summary of the orbital behavior of GX 301-2: hard X-ray flux (a), pulsed fraction (b) and hardness ratio (c). Units along the axis are in solar radii. 



and periodically intercepting the neutron star. The effect of 
the spiral stream accretion on the orbital light curve is the 
pre-periastron flare and the broad bump around the apastron 
phase, both caused by the sudden and large increase in the ac- 
cretion rate. The larger density and narrower geometrical size 
of the intercepted accretion stream near the periastron phase 
causes a flare shorter and more intense than at the apastron 
phase. In this scenario, the same matter that accretes is also 
responsible for a very large absorption column (thus the hard- 
ening) and for scattering of the radiation from the neutron star, 
bringing loss of timing coherence (decrease in the pulsed frac- 
tion). This is indeed what the ASM data show at soft X-rays. 
In the same model, the wind component determines a smooth 
variation of the flux, peaking only at the pre-periastron phase. 
T his is what the SuperAG ILE data show at hard X-rays. 

iLeahv & Kostkal (120081) show that around phase 0.25 the 
soft X-rays are largely absorbed by a partial occultation of the 
neutron star by the atmosphere of the companion star. This is 
consistent with both the ASM flux modulation and with the 
SuperAGILE/ASM hardness ratio. However, we observe an 



even more intense hardening peak of the emitted radiation at 
the pre-periastron that i s not shown in the model prediction by 
ILeahv & Kostkal (12008b . In this case the hardness peak is not 
due to a decrease in the soft X-ray flux, but to a larger increase 
in the hard X-rays than in soft X-rays, still in the presence of 
a strong flare in both energy bands. If absorption is respon- 
sible for such an effect, this implies that the site of the large 
flare lies behind the absorbing material (thus the unabsorbed 
soft X-ray flux would be far larger) and both the flare and ab- 
sorption occur at the same orbital phase, suggesting that the 
two are related in origin. Large increases (factor of ~10) in 
the absorption column ne ar the periastro n phase were indeed 
reported by early studies dHaberllll991l using EXOSAT data, 
lLeahvl[T99ll using TENMA data) while more recent studies, 
using a partial covering absorber mod el, reported smaller vari- 
ations (factor of ~2 with BeppoSA X lLa Barbera et al.ll2005h 
and rather scattered in orbital phase (Mukheriee & Paulll2004 
with RXTE/PCA). 

As already mentioned, at this phase we detect a large de- 
crease in the pulsed fraction, implying that the hard X-rays 
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here are more loosely affected by the NS spinning than the soft 
X-rays. This may suggest that (at least) the hard X-ray emis- 
sion at the pre-periastron peak may be associated with tempo- 
rary conditions setting up when the NS crosses the accretion 
stream at its minimum distance from the companion star. It 
is also interesting to remind that just at this orbital phase, a 
detection of non-thermal radi o emission from GX 30 1-2 has 
been very recently reported by Pes talozzi et al] (120091) . 

Changes of the properties or sites of the hard X-ray emis- 
sion along the orbit are suggested also by our study of the 
neutron star pulse profile as a function of t he orb ital phase. 
Similar to what found by La Bar bera et alj ( 120051) using the 
BeppoSAX data, we detected a large variability in the pulse 
shape of the pulsar in the hard X-ray range. The main points 
of variability we detected are: the structure of the main peak, 
the intensity at the interpulse and the amplitude of the pul- 
sation (i.e., the pulsed fraction discussed above). As found 
by BeppoSAX, in the SuperAGILE data the pulse appears 
smoothest at orbital phases near the pre-periastron. The shape 
we obtained for the orbital phase interval 0.91-0.96 is remi- 
niscent of the BeppoSAX/PDS curve. However, in approach- 
ing and leaving the phase of the large flare, the pulse shape 
shows transient features. A small peak structure (at spin phase 
around 0. 1 in Figure[3) gradually appears on the rising and de- 
caying phases of the flare (orbital phases from 0.91 to 0.06), 
being more prominent after the periastron. Simultaneously, 
the flux level at the inter-pulse phase (around 0.25) becomes 
higher. As the neutron star approaches the apastron and goes 
beyond, the pulse shape at energies above 20 keV becomes 
increasingly complex (although the smaller statistics does not 
allow us to identify individual substructures), with the main 
pulses becoming alternately narrower and wider. Interest- 
ingly, approaching again the pre-periastron phase, a small 
peak appears in the interpulse region. This clearly shows up 
at orbital phase 0.88-0.91. All of these features, including the 
smaller pea k, are not in cont r ast wi th the BeppoSAX results, 
obtained by La B arbera et alj ( 12005b in their detailed analysis 
of the pulse shape in four short orbital phase intervals. Alto- 
gether these results seem to suggest that the emitting region 
on the neutron star is significantly affected by the changing 
accretion pattern along the binary orbit. 

One last consideration derives from the secular variation 
of the pulse period of the neutron star, collected in Figure [4] 
As already com mented and discussed by several authors (e.g. 
iKoh et alJfl997h . the spin history is complex and shows sev- 
eral changes in the spin-up and spin-down trends, sometimes 
also very rapid. The major changes in the past were proba- 
bly accompanied to important variations in the emission. In 
the last ~15 years (from MJD 49500 onward), a general spin- 
down trend is settled, perturbed only by few short and rapid 
spin-up episodes. In particular, despite some small-amplitude 
timing noise that SuperAGILE detected (Figure |4fb)), in the 
last ~two years the GX 301-2 system appears rather steadily 
evolving with a spin-down rate of about 1.3 s yr , showing 
no rapid or large perturbations. 

6. SUMMARY AND CONCLUSIONS 

We studied the long term behavior of the HMXB GX 301-2 
using the hard X-ray data collected by SuperAGILE during 
the first year of the AGILE mission. The observations, car- 
ried out in the 20-60 keV energy band, have a large net ex- 
posure, about 3.7 Ms, and cover all the orbital phases of the 
binary system. The SuperAGILE data offer one of the most 
complete hard X-ray monitoring of the 4 1.5 -day long binary 



period available to date. The data span covers large fractions 
of six orbital cycles. The secular trend of the ~680 s pulse 
period is consistent with the previous observations, although 
the SuperAGILE data show a decrease in t he spin-down rate 
with re spect to the observations reported by lDoroshenko et ail 
(120081) . The spin-down trend is approximately constant, but 
there is indication of a significant timing noise, of the order of 
-0.2%. 

Complementing the SuperAGILE data with those of 
RossiXTEl ASM, the source behavior was characterized at all 
the phases of the binary orbit, in terms of its soft and hard X- 
ray flux, spectral hardness, spin period history, pulsed fraction 
and pulse shape. Compared with the soft X-ray data, the hard 
X-ray lightcurve exhibits a smaller modulation of the emis- 
sion during the orbital phases 0.1-0.8 and a larger intensity 
of the pre-periastron flare. This results in a spectral hardness 
which presents clear orbital modulation, with peaks in corre- 
spondence with the pre-periastron flare and near phase 0.25. 
The timing analysis of the hard X-ray emission showed a vari- 
able pulse shape profile as a function of the orbital phase, with 
substructures detected near the passage at the periastron, and 
a clear modulation of the pulsed fraction, which appears in 
turn strongly anti-correlated with the source intensity. 

Overall, the scenario depicted by the SA data appears gen- 
erally con sistent with the model f or this binary system pro- 
posed by iLeahv & Kostkal ((2008). The predictions of the 
model are based on the variability in the accretion rate caused 
by the orbital motion of the NS in the environment of the mat- 
ter outflow from the companion star. The model is suited to 
the low energy observations by RXTE/ ASM, for which it re- 
produces the peaked flare at the pre-periastron and the broader 
peak at the apastron. They both are related to the NS crossing 
of the accretion stream, the difference between the two being 
mostly in the local density and geometrical size of the stream 
at the crossing point. The higher accretion rate causes a soft 
X-ray brightening, and scattering of the radiation by the same 
matter would cause loss of coherence and then a decrease in 
the pulsed fraction. In the hard X-rays we do not detect the 
broad peak at the apastron phase (although other hard X-ray 
observations with BATSE and SwiftfBAH suggest time vari- 
ability for it). Also, our observation of the folded hardness 
ratio (20-50 keV to 2-12 keV) along the orbit is not fully con- 
sistent with the model expectations. Indeed the model predicts 
an increase in the column density near the apastron phase, but 
it does not predict a similar increase near the pre-periastron, 
that we instead clearly detect. 

More in general, our result s fit in the overall scenario of 
the binary X-ray pulsars (e.g. lParmar et al.ll 19891) . where the 
timing properties are related to the source luminosity. As a 
matter of fact, our analysis shows that the pulsed fraction is 
clearly anti-correlated with the source intensity, allowing us 
to attribute the accretion rate the role of driving parameter 
in the system. However, while the timing properties appear 
rather well correlated with the source intensity, the same is 
not true for the hardness ratio, for which we cannot identify 
a trend with the flux. This implies that the emitting process 
and/or geometry change along the orbit in a way that is not 
immediately associated to the variation of the accretion rate 
only. 
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